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ABSTRACT 

o , 

[ Analyses of the widths and shifts of optically thin emission lines in the ultraviolet 

spectrum of the active dwarf e Eri (K2 V) are presented. The spectra were obtained 
using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and 
the Far Ultraviolet Spectroscopic Explorer. The line widths are used to find the non- 
thermal energy density and its variation with temperature from the chromosphere to 
the upper transition region. The energy fluxes that could be carried by Alfven and 
acoustic waves are investigated, to test their possible roles in coronal heating. Acoustic 
^ waves do not appear to be a viable means of coronal heating. There is, in principle, 

£T} \ ample flux in Alfven waves, but detailed calculations of wave propagation are required 

. before definite conclusions can be drawn about their viability. The high sensitivity and 

spectral resolution of the above instruments have allowed two-component Gaussian fits 
to be made to the profiles of the stronger transition region lines. The broad and narrow 
components which result share some similarities with those observed in the Sun, but 
in e Eri the broad component is redshifted relative to the narrow component and 
contributes more to the total line flux. The possible origins of the two components and 
the energy fluxes implied are discussed. On balance our results support the conclusion 
of Wood, Linsky & Ayres, that the narrow component is related to Alfven waves 
reaching to the corona, but the origin of the broad component is not clear. 

Key words: stars: individual (e Eridani) - stars: late-type - stars: chromospheres - 



stars: coronae - line: profiles 



1 INTRODUCTION 

The observed widths (and profiles) of emission lines, formed 
in stellar chromospheres, transition regions and coronae, 
provide valuable information on possible heating processes, 
in addition to that deduced from the emission line fluxes. 
Observations of solar line widths in the ultraviolet (uv) part 
of the spectrum have been carried out since 1970 (Berger, 
Bruner & Stevens 1970; Bruner et al. 1970) and much of the 
basic behaviour of these line widths was established from 
work in the 1970s using spectra obtained from rocket flights 
(Boland et al. 1973, 1975; Kjeldseth-Moe & Nicolas 1977) 
and from the Naval Research Laboratory's (NRL) spectro- 
graph on Skylab (Doschek et al. 1976; Feldman et al. 1976). 
Further work was carried out with instruments on Orbiting 
Solar Observatory IV (OSO-IV) and the Solar Maximum 
Mission (SMM). These early observations showed that line 
broadening is present in excess of the thermal Doppler width 
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and that when interpreted in terms of a wave energy flux, 
this flux is essentially constant through the transition re- 
gion from T e ~ 10 4 K to ~ 10 5 K (where T e is the electron 
temperature). The observed line widths were found to be 
very similar in the various solar locations studied. Using 
spectra obtained with the NRL's High Resolution Telescope 
and Spectrograph (HRTS), Kjeldseth-Moe & Nicolas (1977) 
found that some emission lines could not be well-fitted with 
a single Gaussian profile, but have excess flux in a broad 
component. Such profiles were studied further by Dere & 
Mason (1993). The earliest instruments used photographic 
recording, which has some limitations, but the main results 
have been confirmed and extended using the Solar Ultravi- 
olet Measurements of Emitted Radiation (SUMER) instru- 
ment on the Solar and Heliospheric Observatory (SOHO). 
This has allowed line profiles to be studied with a combina- 
tion of high sensitivity and high spatial and spectral resolu- 
tion, at a variety of locations (see e.g. Doyle et al. 1997; Seely 
et al. 1997; Chae, Schiihle & Lemaire 1998). Peter (2000a,b; 
2001) has made systematic measurements of line widths and 
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shifts in both supergranulation cell boundary and cell inte- 
rior regions. We make comparisons with the results of Chae 
et al. (1998) and Peter (2000a,b; 2001) in later sections of 
this paper. 

Similar studies of line widths and line shifts were made 
for a few late-type stars observable at high resolution with 
the International Ultraviolet Explorer (WE) (see e.g. Ayres 
et al. 1983a,b; Brown et al. 1984; Jordan et al. 1987). Signifi- 
cant improvements over these were obtained using the God- 
dard High Resolution Spectrograph (GHRS) on the Hub- 
ble Space Telescope (HST), and, most recently, using the 
greater spectral coverage and sensitivity of the Space Tele- 
scope Imaging Spectrograph (STIS). The observations with 
IUE showed the presence of line broadening in G/K-type 
dwarfs which is similar to that in the Sun, and broader line 
widths in evolved stars with hot coronae. Observations of 
dwarf and giant stars with the GHRS have shown the pres- 
ence of excess flux in line wings, compared with single Gaus- 
sian fits to the profiles (Linsky & Wood 1994; Wood, Linsky 
& Ayres 1997, and references therein); the contribution of 
this broad component appears to increase with indicators of 
stellar activity. 

The causes of the observed line broadening are not 
known, but several possible processes have been proposed. 
The excess width of the lines, compared with the thermal 
broadening, is commonly attributed to turbulence, and ob- 
served line profiles have been used to find the most-probable 
turbulent velocities or the root-mean-square (r.m.s.) veloci- 
ties. These have been used to find non-thermal energy densi- 
ties and the energy flux carried by the non-thermal motions, 
assuming that they are associated with the passage of waves 
through the atmosphere. (See e.g. Brown et al. 1984; Jor- 
dan et al. 1987; Wood et al. 1997 for stellar applications). 
When non-Gaussian profiles have been observed, the excess 
flux in the wings has been attributed to a broad component 
associated with 'explosive events' ('microflaring') (Linsky & 
Wood 1994), or, in the Sun, to the passage of Alfven waves 
(Peter 2001). The effect of acoustic waves on line profiles 
has also been discussed (Bruner & McWhirter 1979). 

The possible origins of waves in the solar atmosphere 
and the processes by which these and other mechanisms, 
in particular magnetic reconnection, may cause heating of 
the solar chromosphere and corona have been the topic of 
many papers and reviews. We have made use of reviews by 
Browning (1991), Narain & Ulmschneider (1990, 1996) and 
Zirker (1993). 

Here we report observations of emission line profiles, 
widths and line shifts in the active dwarf e Eri (K2 V) . These 
have been made with the STIS and with the instruments on 
the Far Ultraviolet Spectroscopic Explorer (FUSE). The data 
presented here are the most complete so far published for a 
main-sequence star, other than the Sun. 

The new observations and data reduction are described 
in Section 2. Analyses in terms of single Gaussian fits to the 
line profiles and the implications of the line widths for the 
energy fluxes carried by Alfven waves and acoustic waves 
are discussed in Section 3. The results of the measurements 
of line shifts are given in Section 4. Some of the line profiles 
are fitted better by 2-component profiles rather than single 
Gaussian fits; the widths and shifts derived for the two com- 
ponents and their possible interpretations are described in 
Section 5. Our conclusions are summarized in Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 

The STIS observations of e Eri have been described by Jor- 
dan et al. (2001a). The observations made with FUSE were 
obtained from the archive. In using the latter, we restrict 
our analyses to the lines of Cll and Ovi, since these were 
obtained with the la lithium-fluoride detector, which is the 
one used to obtain the pointing on the star. The emission 
lines studied and the transitions involved are listed in Table 
1, together with their rest wavelengths (Ao), observed wave- 
lengths (A obs ) and line widths (AA), which are the full- width 
at half- maximum (FWHM) . The rest wavelengths are taken 
from Kurucz and Bell (1995) or from Griesmann & Kling 
(2000), with the exception of that for the Fe XII line, which 
is based on the solar wavelength (see Jordan et al. 2001a). 
The observed wavelengths and comparisons with the rest 
values given by the above authors and by Kelly (1987) are 
discussed in Section 4. 

The widths given in Table 1 were measured by mak- 
ing single- Gaussian fits to the observed spectra. Except at 
wavelengths above ~ 1800 A, the background level is close to 
zero. At A > 1800 A, the background was determined from 
the nearby continuum. Likely errors (based on the wave- 
length resolution of the appropriate echelle grating, the un- 
certainty in the background level, the goodness of the fit 
and comparison with fluxes found by direct integration) are 
given for each measurement. For some of the weakest lines 
the signal-to-noise ratio was judged to be too small to give 
reliable line widths. In addition to the Gaussian fits, for the 
stronger lines, interactive measurements of the width at half 
the peak amplitude were made. For most of the lines these 
agree with the Gaussian fit values, within the error bars, but 
for the resonance lines of Civ, Nv, Ovi and Siiv they are 
significantly smaller indicating that these profiles are not 
adequately described by a single Gaussian fit. These half 
peak-amplitude widths are also given in Table 1. It is found 
that the profiles of these lines are described much better by 
two-Gaussian fits (see Section 5). 

The STIS instrumental widths are taken from the STIS 
WWW site, and are 1.1 - 2 pixels, depending on the wave- 
length and grating used. For the E140M grating, one pixel is 
typically ~ 0.015 A, and this instrumental broadening does 
not contribute significantly to the observed line widths. For 
the E230M grating, one pixel is typically ~ 0.03 A, which 
leads to significant contributions to the widths of longer 
wavelength lines observed with this grating (see Table 1 in 
Jordan et al. 2001b). The FUSE observations have a resolu- 
tion R = 20,000 (Blair & Andersson 2001), corresponding 
to AA ~ 0.09 A. 

The effect of rotational broadening has been considered 
but is judged to be negligible. For e Eri, vsini < 2 km s _1 
(Saar & Osten 1997) which is too small (compared to the in- 
strumental resolution) to be detectable; it is noted however, 
that only a factor of a few increase in the spectral resolution 
would make this detectable. 



3 NON-THERMAL VELOCITIES 

Most of the permitted lines, and all of the intersystem lines, 
formed in the transition region are expected to be optically 
thin. In the transition region, line opacities can be estimated 
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Table 1. Transitions, wavelengths and widths of the emission lines studied. 



Ion 



Transition 



Ao (A) A obs (A) 



AApwHM (A) 



C i 
C i 
C i 
C in 
C in 
C in 
C in 
C iv 

C iv 

N iv 
N v 

N v 

O i 
O i 
O in 
O in 
O iv 
O iv 
O iv 



2s 2 2p 2 3 P 



S i 
S i 
S ii 
S ii 
CI i 
Fe ii 
Fe ii 
Fe ii 
Fe II 
Fe II 
Fe II 
Fe XII 



- 2s2p 3 3 Di 
2p 2 3 Pi - 2p3s 3 P 2 
2p 2 X D 2 - 2p3s 3 Pi 
2s2p 3 Pi - 2p 2 3 P 2 
2s2p 3 P - 2p 2 3 Pi 
2s2p 3 Pi - 2p 2 3 P 



2s2p 3 P 2 
2s 2 S 1/2 - 



2 3 



2p 

2p 2 P 



Pi 

3/2 



2s 2 S 1/2 -2p 2 P 1/2 

2s 2 1 S - 2s2p 3 Pi 
2s 2 S 1/2 -2p 2 P 3/2 

2s 2 S 1/2 - 2p 2 P 1/2 

2p 4 3 P 2 -2p 3 ( 4 S)3s 5 S 2 
2p 4 3 Pi -2p 3 ( 4 S)3s 5 S 2 
2s 2 2p 2 3 Pi - 2s2p 3 5 S 2 
2s 2 2p 2 3 P 2 - 2s2p 3 5 S 2 
2s 2 2p 2 P 1/2 - 2s2p 2 4 P 1/2 



2s 2 2p 2 P ; 



2s 2 2p 2 P. 



3/2 
3/2 



2s2p 2 4 P 
2s2p 2 4 P 



/2 
1/2 



O v 


2s2p J Pi 


- 2p 2 *D 2 


O vi 


2s 2 S 1/2 


- 2p 2 P 3/2 


O vi 


2s 2 S 1/2 


- 2p 2 P 1/2 


Al ii 


3s 2 % - 


3s3p !Pi 


Si ii 


3s 2 3p 2 P 3/2 


- 3s3p 2 2 P 3/2 


Si ii 


3s 2 3p 2 P 3/2 


- 3s3p 2 2 P 1/2 


Si ii 


3p 2 P 3 /2 


- 3d 2 D 5/2 


Si ii 


3p 2 P 3 / 2 


- 3d 2 D 3/2 


Si ii 


3p 2 Pi/ 2 


-4s 2 S 1/2 


Si ii 


3p 2 P 3/2 


-4s 2 S 1/2 


Si ii 


3s 2 3p 2 P 1/2 


- 3s3p 2 2 D 3/2 


Si ii 


3s 2 3p 2 P 3/2 


- 3s3p 2 2 D 5/2 


Si ii 


3s 2 3p 2 P 3/2 


- 3s3p 2 2 D 3/2 


Si ii 


3s 2 3p 2 P 1/2 


- 3s3p 2 4 P 1/2 


Si ii 


3s 2 3p 2 P 3/2 


- 3s3p 2 4 P 1/2 


Si in 


3s3p 3 Pi 


- 3p 2 3 P 2 


Si in 


3s3p 3 P 


- 3p 2 3 Pi 


Si in 


3s3p 3 Pi 


- 3p 2 3 P 


Si in 


3s3p 3 P 2 


- 3p 2 3 Pi 


Si in 


3s 2 l S - 


3s3p 3 Pi 


Si iv 


3s 2 S 1/2 - 


- 3p 2 P 3/2 


Si iv 


3s 2 S 1/2 - 


- 3p 2 P 1/2 



3p 4 3 P 2 
3P 4 " 



3p 3 ( 4 S)4s 5 S 2 
3p 3 ( 4 S)4s 5 S 2 



3s 2 3p 3 4 S 3/2 -3s3p 4 4 P 3/2 
3s 2 3p 3 4 S 3/2 -3s3p 4 4 P 5/2 
3p 5 2 Pi/ 2 - 3p 4 ( 3 P)4s 2 P 1/2 
3d 6 (a 5 D)4s a 6 D 9/2 - 3d 5 (a 6 S)4s4p( 3 P) z 8 P 9/2 
3d 6 (a 5 D)4s a 6 D 9/2 - 3d 5 (a 6 S)4s4p( 3 P) z 8 P 7/2 
3d 6 (a 5 D)4s a 6 D 7/2 - 3d 5 (a 6 S)4s4p( 3 P) z s P 7/2 
3d 6 (a 5 D)4s a 6 D 7/2 - 3d 5 (a 6 S)4s4p( 3 P) z 8 P 5/2 
3d 6 (a 5 D)4s a 6 D 5/2 - 3d 5 (a 6 S)4s4p( 3 P) z 8 P 5/2 
3d 6 (a 5 D)4s a 6 D 3/2 - 3d 5 (a 6 S)4s4p( 3 P) z 8 P 5/2 
3p 3 4 S 3/2 - 3p 3 2 P 3/2 



1560.310 
1656.267 
1993.620 
1174.933 
1175.263 
1175.987 
1176.370 
1548.204 



1486.496 
1238.821 



1355.598 
1358.512 
1660.809 
1666.150 
1399.780 
1401.157 
1407.382 
1371.296 
1031.912 



1670.789 
1194.500 
1197.394 
1264.738 
1265.002 
1526.707 
1533.431 
1808.013 
1816.928 
1817.451 
2335.123 
2350.891 
1294.545 
1296.726 
1301.149 
1303.323 
1892.030 
1393.760 



1900.287 
1914.697 
1253.811 
1259.519 
1351.656 
1888.010 
1901.773 
1915.792 
1926.240 
1936.794 
1944.134 
1242.00 



1560 
1656. 
1993 
1174 
1175 
1176 
1176 
1548 



331 ± 0.008 
273 ± 0.009 
642 ± 0.016 
949 ± 0.006 
281 ± 0.006 
002 ± 0.006 
384 ± 0.006 
224 ± 0.008 



1550.781 1550.788 ± 0.008 



1486.521 ± 0.017 
1238.828 ± 0.007 



1242.804 1242.822 ± 0.007 



1355.607 
1358.523 
1660.847 
1666.169 
1399.797 
1401.186 
1407.429 
1371.327 
sec text 



± 0.007 
± 0.007 
± 0.018 
± 0.018 
± 0.015 
± 0.008 
± 0.015 
± 0.007 



1037.613 see text 



1670. 
1194. 
1197. 
1264. 
1265. 
1526. 
1533. 
1808. 
1816. 
1817. 
2335. 
2350. 
1294. 
1296. 
1301. 
1303. 
1892. 
1393. 



802 ± 0.009 
499 ± 0.013 
389 ± 0.013 
747 ± 0.007 
012 ± 0.007 
717 ± 0.008 
435 ± 0.008 
069 ± 0.015 
970 ± 0.015 
491 ± 0.015 
109 ± 0.019 
896 ± 0.019 
578 ± 0.014 
758 ± 0.014 
168 ± 0.028 
349 ± 0.014 
106 ± 0.015 
781 ± 0.008 



1402.773 1402.795 ± 0.008 



1900. 
1914. 
1253. 
1259. 
1351. 
1888. 
1901. 
1915. 
1926. 
1936. 
1944. 
1241. 



322 ± 0.016 
757 ± 0.016 
820 ± 0.015 
533 ± 0.015 
664 ± 0.007 
032 ± 0.015 
799 ± 0.016 
817 ±0.016 
278 ± 0.016 
810 ± 0.016 
156 ± 0.016 
992 ± 0.007 



0.135 ±0.008 6 
0.191 ±0.008 6 
0.110 ±0.016 
0.157 ±0.015 
0.155 ±0.015 



0.207 ± 
0.183 ± 
0.210 ± 
0.177 ± 
0.210 ± 
0.178 ± 
0.155 ± 
0.172 ± 
0.149 ± 
0.065 ± 
0.056 ± 



0.008 

0.008° 

0.008 

0.008° 

0.050 

0.007 

0.007 a 

0.007 

0.007 a 

0.007 

0.010 



0.217 ±0.030 



0.191 ±0.015 



0.193 ±0 
0.194 ±0 
0.180 ±0. 
0.185 ±0 
0.165 ±0 
0.213 ±0. 
0.188 ±0 

_C 

0.184 ±0. 
0.157 ±0 
0.173 ±0 
0.181 ±0. 
0.136 ±0. 
0.142 ±0 
0.117±0 
0.129 ±0. 
0.124 ± 



0.194 ± 
0.176 ±0. 
0.143 ±0 
0.175 ±0. 
0.155 ±0 
0.104 ±0 
0.106 ±0 
0.108 ±0. 
0.122 ±0. 
0.061 ±0 



011 

015 

015 a 

015 

015 a 

ois 6 '' 

040 6 

007 6 

007 6 

008 6 

008 6 

015 b ' ! 

015 6 '' 

015^ 

019 

019 



015 

008 

008 a 

008 

008 a 

016 

016 

015 

015 

007 



0.130 ±0.016 

_C 

0.127 ±0.016 
0.124 ±0.016 

_C 

0.21 ±0.04 



a Width at half the peak amplitude. Optically thick line. c Line too weak to determine the width. 
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from the observed line fluxes and the mean electron pres- 
sure determined by Jordan et al. (2001b). Intersystem lines 
formed in the chromosphere are also estimated to be opti- 
cally thin. It is harder to estimate the opacities of permitted 
chromospheric lines and, as will be discussed later, some of 
these are likely to be optically thick, and will have contribu- 
tions to their line widths from photon scattering. Lines to 
the ground level of atoms and first ions may also be affected 
by interstellar absorption. The ratios of lines within multi- 
plets and the presence of self-reversals have also been used 
to judge whether or not lines are optically thin. 

As discussed in Section 1, the widths of optically thin 
emission lines in stellar uv spectra are significantly greater 
than the thermal Doppler width corresponding to the tem- 
perature of line formation (Tf orm ). When corrected for in- 
strumental broadening it is usually assumed that the ob- 
served width is the result of thermal and non-thermal mo- 
tions, the latter being referred to as turbulence. Since most 
optically thin line profiles are well described by one (or more) 
Gaussian components, it is assumed that the turbulent ve- 
locities follow a Maxwellian distribution (and hence produce 
a Gaussian line profile) , and the distribution is characterised 
by its most probable velocity £. This non-thermal contribu- 
tion to the line width is usually the dominant contribution 
to chromospheric and transition region line widths. 

The total line width may be expressed as 



Table 2. Values of Tf orm adopted. 



V Ao / obs 



/AA\ 2 (-Y (—\ 2 

V Ao J 

therm V Ao / inst V Ao / 

non-them 



(1) 



where the subscript "obs" , "therm" , "inst" and "non-therm" 
refer to the observed total width and the contributions from 
thermal, instrumental and non-thermal broadening respec- 
tively. By expressing the thermal width in terms of the tem- 
perature of line formation and the non-thermal width in 
terms of a most probably velocity £, equation (1) may be 
written as 



41n2 



giving 

e 



= 0.361 



/AA\ 2 _ / AA\ 2 

\ Ao / obs V Ao / inst 



/AA\ 2 _ / AA\ 2 

V Ao / obs V Ao / inst 



2/CB?form 

m H mionC 2 



(2) 



1.84xl0- 13 ^^(3) 



where Tform is in Kelvin, mu is the mass of hydrogen and 
mion is the mass of the emitting ion, relative to that of hy- 
drogen. 

For the chromospheric lines, formed at low tempera- 
tures, the value of £ is relatively insensitive to the adopted 
value of Tform, but because these lines are formed over a 
wide range of temperatures, the actual values of Tfo rm arc 
less meaningful. For the transition region lines the value of £ 
is sensitive to Tf or m, which should be regarded as only an av- 
erage value. The adopted values of Tform are given in Table 2. 
For the transition region lines formed above log T e > 4.3, the 
values of Tform are based on the line contribution functions 
and the observed emission measure distribution (Jordan, 
Sim & McMurry 2001c; Sim 2002). New calculations of the 
relative ion populations by Sim (2002) have been used, which 
take into account the density sensitivity of di-electronic re- 
combination. This reduces the temperature of optimum line 
formation compared with the zero density calculations of Ar- 
naud & Rothenflug (1985). Overall, the values of Tform arc 



Ion 


Ao (A) 


log T form (K) 


C I 


Ty yo . 


3.7 


C in 


1 1 70 mult 

± L 1 U IllUllr. 


4.7 


C iv 


1 ^fl mult 

llJUU IIlLlllr. 


4.9 


N iv 




5.0 


N v 


1 940 mult 


5.2 


O i 


1356 miiltr. 


O.o 


O in 




4.9 


O iv 


1401.16 


5.2 


O v 


1371.30 


5.4 


O vi 


1034 mult. 


5.45 


Al ii 


1670.79 


4.2 


Si ii 


1815 mult. 


3.85 


Si ii 


2340 mult. 


3.65 


Si in 


1892.03 


4.45 


Si iv 


1400 mult. 


4.8 


S i 


1907 mult. 


3.7 


S n 


1255 mult. 


4.2 


CI I 


1351.66 


3.7 


Fe II 


1915 mult. 


3.85 


Fe xii 


1242.00 


6.15 



expected to be accurate to within about ±0.1 dcx. Values 
for the chromospheric lines of O I, C I and Si n are based on 
Sim's (2002) calculations with the MULTI (Carlsson 1986, 
Scharmer & Carlsson 1985) radiative transfer code; those 
for CI I, Al n and Fe II are only rough estimates; the value 
adopted for Fe xii is that used by Jordan et al. (2001a). 

The values of £ found using equation (3) are given in Ta- 
ble 3. The errors on £ include both the measurement errors 
from Table 1 and an uncertainty in logTfo rm of ±0.1 dex. 
The value in parenthesis after the total error is the contribu- 
tion from the uncertainty in log Tfo rm alone. The observed 
line positions have been converted to velocity shifts using 
the rest wavelengths Ao given in Table 1. 

The values of £ given in Table 3 are plotted against 
Tform in Fig. 1, including only those lines which are esti- 
mated to be optically thin. Qualitatively, this plot is similar 
to those found for the Sun (see e.g. Jordan 1991; Chae et 
al. 1998); £ is smallest in the chromosphere, rises in the low 
transition region and then flattens off. Comparing optically 
thin lines in common with those of Chae et al. (1998), the 
range of £ is slightly smaller in e Eri, the widths of the mid- 
transition region lines are smaller in e Eri, but the widths 
of the neutral lines are on average larger in e Eri. All the 
lowest temperature chromospheric lines (those of the neu- 
trals) are consistent with £ ~ 7.5 km s _1 . The widths of the 
lowest temperature lines of Si II (the intersystem multiplet 
around 2340 A) a gree with this mean value. All the lines with 
logTform > 4.5 have widths that are consistent with a con- 
stant value of £ m 21.3 km s _1 throughout the mid-transition 
region. The other lines indicate a smooth transition between 
the two regimes. 

There are some differences in the value of £ for lines 
formed at similar temperatures, which may be due to opac- 
ity effects. The value of £ for the Al II resonance line is larger 
than those for the S II lines. It is difficult to judge whether or 
not the observed Al II profile is self-reversed, or contains in- 
terstellar absorption, but the opacity estimate suggests that 
small optical depth effects could occur. The lines in the Si II 
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Table 3. Derived line shifts (u) and turbulent velocities (£). The error estimates for £ include the total uncertainty due to the measurement 
errors and an assumed uncertainty in logTf orm of ±0.1 dex. The errors in parentheses indicate the uncertainty due to AlogTf orm alone. 



Ion 


Ao (A) 


v (km s ) 


£; (Km S ) 


Ion 




ii -l\ 
v (km s ) 


i; (Kill S ) 


C i 


1560.310 


+4.0 ± 1.5 




Si ii 


1526.707 


+2.0+1.6 


— 


C i 


1656.267 


+1.1 ± 1.6 




Si ii 


1533.431 


+0.8 + 1.6 




C i 


1993.620 


+3.3 + 2.4 


+ 1 90C0 101 
'•° -2.20(0.12) 


Si ii 


1808.013 


+9.3 + 2.5 


„ „ 7(1(0 041 
iZ * U -1.77(0.05) 


C in 


1174.933 


+4.1 ± 1.5 


_ „ _„_i_2 Qd(0 45) 
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multiplet around 1815 A do not have ratios which agree with 
their optically thin values, and the two lines with the largest 
oscillator strengths (1808.01 A and 1816.93 A) are slightly 
broader than the line at 1817.45 A and the intersystem lines 
of Fe ii. Although these Si n lines are not self-reversed, small 
opacity effects could be occurring. For these reasons these 
lines of Si n and Al n are not used in subsequent analyses 
using non-thermal widths. Lines from the other multiplets of 
Si n have significantly larger oscillator strengths than those 
in the 1815-A multiplet, and should therefore show opacity 
effects. The lines at 1533.43 A and 1526.71 A have a flux 
ratio which is ~ 1, instead of the optically thin value of 2. 
The multiplets at shorter wavelengths have larger or simi- 
lar oscillator strengths (as given in CHIANT1, Landi et al. 
1999), and are therefore also expected to be optically thick. 
The widths of these other lines of Si n are included in Table 
1, but are not included in Fig. 1. 



3.1 Kinetic energy density 

The interpretation of the line broadening in spatially inte- 
grated stellar spectra is subject to some uncertainties. Solar 
observations can be used to investigate the behaviour of line 
widths at different positions on the Sun. In both early stud- 
ies (see references in Section 1) and in more recent work 
with SUMER (see e.g. Peter 1999), when single Gaussian 
fits are made to the line profiles there is little systematic dif- 
ference between the widths at different solar locations and 
distances from the limb. In this respect, using spatially av- 
eraged spectra should give a good estimate of the typical 
line widths. These solar observations also show that there 
are no preferred directions of the motions contributing to 
the widths. However, in the Sun, there are line (Doppler) 
shifts, with (on average) red-shifts being observed in the su- 
pergranulation network regions and smaller red-shifts being 
observed in cell interior regions (Peter, 2000a, b). (Blue-shifts 
are also observed at individual locations- see Section 4.) In- 
tegrated over a star, the combined effects of these motions 
will contribute to the observed line widths. Peter (2000a,b) 
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Figure 1. log£ versus logTf orm . The vertical error bars indicate 
the total uncertainty in the deduced values of £ (including mea- 
surement errors and uncertainty in Tf orm ) and the sloping bars 
indicate the effects of uncertainties in Tf orm . 



Figure 2. The non-thermal energy density in the gas plotted as 
log {/nt versus log Tf orm . The vertical bars show the uncertainty 
in C/nt due to the total uncertainty in £, but make no allowance 
for possible errors in p. 



has made quantitative estimates of this effect by studying 
widths and shifts from relatively small regions with those in 
the mean spectrum observed over a much larger area. Using 
Peter's (2000b) corrected widths, the values of £ deduced 
from the mean spectra would overestimate the median val- 
ues by about 9 per cent (narrow component) to 15 per cent 
(broad components). (Two-component fits to our spectra are 
discussed in Section 5). For all except the stronger (alkali- 
like) lines, these potential contributions are less than our 
estimated uncertainties in the line widths. 

Assuming that the line widths are due to non-thermal 
motions they can be used to estimate the total non-thermal 
kinetic energy density in the gas, L r NT(7 1 e ), which is given 
by 



E/nt = ^p£ 2 



(4) 



where p is the mass-density. This follows from the as- 
sumption that the non-thermal velocities follow a Maxwell- 
Boltzmann distribution in three dimensions, with most 
probably speed £. (See the Appendix for the cases of mo- 
tions in one and two dimensions), p is given by 



p = m,Hp.Ng 



(5) 



where ran is the mass of hydrogen, N g is the number density 
of all particles and p is the mean atomic weight defined as 



I" = 



J2 i miNi 



(6) 



where rrii and Ni are the mass and number density of par- 
ticles of species i in the plasma and the sums run over all 
such species. Adopting abundances similar to those in the 
solar photosphere, p can be approximated by 



/' : 



1.4/Vh 



iV c + ITTVh 



(7) 



where Nn is the number density of hydrogen and N e is the 
electron number density. This gives pN g « 1.4/Vh. Hence 



[7nt « l.lra H iVH£ 



(8) 



Fig. 2 shows a plot of [/nt in the gas against temper- 
ature. Nh has been taken from a semi-empirical model of 
the average atmosphere of e Eri (Sim 2002). The model is 
based on an emission measure distribution derived from line 
fluxes observed with STIS and the chromospheric model 
of Thatcher, Robinson & Rees (1991). This new model, 
and others, will be discussed in a forthcoming paper (Sim, 
in preparation). The important parameters in the model 
adopted here (column mass density, electron temperature, 
electron pressure, proton pressure, hydrogen pressure and 
sound speed) are shown in Fig. 3. The error bars in Fig. 2 
give the total uncertainty in the values of £, but not the un- 
certainties in p. The form of C/nt with T e is dominated by 
the variation of p, and is therefore considerably more sensi- 
tive to the adopted atmospheric model than to the measured 
values of £ (p varies by more than five orders of magnitude 
across the temperature range in Fig. 2, compared to only 
a factor of three variation in £). It is clear, however, that 
the non-thermal energy density in the chromosphere is at 
least three orders of magnitude greater than in the overly- 
ing transition region and corona. 

When considering the propagation of waves which are 
coupled to the magnetic field (torsional Alfven waves and 
fast mode waves when the Alfven speed greatly exceeds the 
sound speed), the degree of ionization of the gas needs to 
be considered. This is relevant to the chromospheric regions 
where hydrogen is not fully ionized (see e.g. discussion in 
Osterbrock 1961). If the frequency of collisions between ions 
and neutrals is much larger than the wave frequency then 
there will be collisional coupling between the ions and neu- 
tral atoms, which can cause damping of the waves. In this 
case the gas density p is appropriate. Conversely, if the col- 
lision frequency is much less than the wave frequency, the 
ion density p ion becomes the relevant density. The chromo- 
spheric model has been used to calculate pion, allowing for 
the protons, 'metals' with low first ionization potential, and 
the electrons produced by these and hydrogen (i.e. the ions 
and electrons are treated as a single fluid) . p ion is then given 
by 
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Figure 3. The adopted atmospheric model. The left-hand panel shows the electron temperature (T e ) versus the mass column density 
(ficol)- The right-hand panel shows the electron (solid line), proton (dotted line) and total hydrogen (H and H + , dashed line) pressures, 
and also the sound speed (v B , triple-dot dashed line) versus T c . 
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Figure 4. The non-thermal energy density in the ions plotted as 
log [/fjx versus log Tf orm . The vertical bars show the uncertainty 
in C/nt due to the total uncertainty in £, but make no allowance 
for possible errors in Pi on - 



m c N e + mpN p + rn m (N e — N p ) 



(9) 



where m c and m p are the masses of the electron and proton, 
respectively, m m is the mean mass of ions of elements other 
than hydrogen divided by the mean ionisation stage of such 
ions, and N p is the proton number density. p ion then replaces 
p in equation (4) . 

Fig. 4 shows the non-thermal energy density in the ions, 
which at T e < 10 4 K is significantly smaller than the total 
non-thermal energy density shown in Fig. 2. Since collisions 
between the ions and neutrals will occur and will transfer 
some energy to the neutrals, the actual relevant non-thermal 
energy density will lie between the extremes given by Figs. 
2 and 4. 



3.2 Wave Heating 

We now consider the possibility that the non-thermal energy 
density is associated with waves propagating at either the 



Alfven speed or the sound speed, and examine the extent to 
which these hypotheses can account for the observations. 

Provided the wave frequencies are sufficiently small (so 
that the WKB approximation is valid), and non-linear ef- 
fects are small, the non-thermal energy flux, Fnt, is related 
to the kinetic energy density by 



(10) 



where v prop is the appropriate propagation velocity and the 
factor of 2 accounts for equipartition of energy between ki- 
netic and potential energies in the waves. 

First we consider waves which propagate at the Alfven 
speed, without distinguishing between torsional (shear) 
waves or compressional (fast-mode) waves when v a >> v s 
(see Cross 1988 for Alfven wave theory and references to 
early work). The formulation is given in terms of Pi on - 



B 



(11) 



where B is the magnetic flux density. B is not known directly 
above the photosphere. Riiedi et al. (1997) have measured 
the spatially averaged surface magnetic flux density, Bo = 
B S A S /A t , where B S A S is the surface magnetic flux. They 
find a value of 165 ± 30 G. In the overlying atmosphere we 
make the assumption that some fraction of the total stellar 
surface area is occupied by flux tubes, each with the same 
values of B and individual cross-sectional areas. Since not 
all the surface magnetic flux may continue through to higher 
layers, we specify the fraction which continue at a given T e 
as b = B A/ Bo At,, where A* is the total surface area and A is 
the actual total surface area occupied by the magnetic field 
B, at given T e . The observed spatially averaged non-thermal 
flux is then 



p A 3 2 A 

~A = ^Piont, Va-T- 



2V47T 



l/2 f 2. D -2 -1 

P ion i bB o erg cm s 



(12) 



where F NT refers to an individual flux tube. The factor 3, 
which was introduced in the calculation of the non-thermal 
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Figure 5. log F^tA/A* versus log Tf orm for Alfven waves, adopt- 
ing 6=1. The vertical error bars indicate the uncertainty in -Fnt 
due to the total uncertainty in £; no allowance is made for uncer- 
tainties in the values of pi on in the atmosphere. The open symbols 
indicate the calculations using pi on and the filled symbols show 
those with the total density p. The full line shows the flux required 
to account for the radiation losses above a given T e . 



energy density by the assumption of isotropic non-thermal 
motions in three dimensions, has been retained here. If the 
motion were in fact restricted to only 2 or 1 dimensions, 
as associated with purely transverse or longitudinal oscilla- 
tions, equation (12) would overestimate the actual flux by 
factors of 1.5 and 3 respectively. 

The spatially averaged non-thermal flux calculated us- 
ing equation (12) and the observed value of Bo = 165 G 
is shown in Fig. 5, initially adopting 6=1. Results using 
p and pi on are illustrated. The flux required to account for 
the radiation losses above a given T e is also shown. The ra- 
diation losses have been calculated using the model shown 
in Fig. 3, which is based on the spatially averaged emission 
measure distribution derived by Sim (2002) and the radia- 
tive power losses of Cook et al. (1989), adopting the stellar 
photospheric abundances (which are similar to those of the 
Sun). 

We also consider heating by acoustic waves propagating 
at the sound speed, given by 



-yk B T e 
pmH 



1/2 



(13) 



where 7 = 5/3 is the ratio of specific heats and T e is the 
electron temperature. The mean molecular weight is given 
by equation (7) . If the acoustic waves originate in some frac- 
tion of the surface area (A /A*), the corresponding spatially 
averaged non-thermal energy flux is 



A_ 



A 



3 t a 
2* Vs A 



(14) 



4.1 x 10 



p A, 



erg cm 



Fig. 6 shows the non-thermal flux assuming that the 
turbulent widths are due to the passage of acoustic waves 
with A /A* = 1, and the flux required to account for the 
radiation losses above a given T e . Above 10 4 K, the forms of 



Figure 6. log Fnt versus log Tf orm for acoustic waves. The ver- 
tical error bars indicate the uncertainty in Fnt due to the total 
uncertainty in f, but do not attempt to estimate the uncertainty 
in the variation of p in the atmosphere. The full line shows the 
flux required to account for the radiation losses above a given T e . 



Figs. 5 and 6 are very similar; this is because, apart from the 
numerical constant, the only difference betwe en equations 
(12) and (14) is the factor y/T e (N e + I.INh) = yfP g . The 
pressure scale height is much greater than the thickness of 
the transition region and so P g varies little between log T e = 
4.0 and logT e = 6.0. However P g increases rapidly below 
T e ~ 10 4 K, which is why Fig. 6 displays a steeper gradient 
in the chromosphere. 

The form of the curves in Figs. 5 and 6 show that in 
both cases the non-thermal wave flux is roughly constant 
between logTe = 4.2 and 4.8, consistent with a transport 
of non-thermal energy with little dissipation. This apparent 
constancy of the wave flux also shows up in Fig. 1, since 
over this range of temperature £ is approximately tx T^ A , 
as in the solar transition region (Jordan 1991). This gradient 
is expected if there is a constant Alfven wave flux (in the 
WKB approximation) , or a constant acoustic flux, provided 
the electron pressure and b are constant. Under the same 
circumstances a constant turbulent energy flux would give 
£ oc T^ 3 . At higher temperatures, £ is roughly constant 
and then decreases, so Fnt appears to decrease. This could 
occur if energy is being dissipated (but see below). It is 
not surprising that £ hardly varies over the mid-transition 
region. The temperature gradient in the spatially averaged 
models is large (as in the Sun), so the range of T e involved 
corresponds to a small height difference (less than 100 km) . 

In the case of an acoustic energy flux it is possible that 
the decrease above 10 5 K is simply a geometrical effect; the 
area occupied by the regions associated with the acoustic 
flux could increase with T e . The flux shown in Fig. 6 is the 
spatially averaged value, with the term Aj 'A* on the r.h.s of 
equation (14) taken as 1. An increase in A /A* from a low 
value in the transition region to a larger value in the corona 
could lead to FntA/A, being constant until dissipation oc- 
curs in the corona. 

Simple energy requirement arguments can probably 
rule out coronal heating by acoustic waves. From Fig. 6 
the acoustic energy flux passing into the corona at T e ~ 
1.4 x 10 6 K is w 8 x 10 5 ergs cm~ 2 s _1 , assuming that 
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Figure 7. log£ versus logTf orm from the solar data of Chae et 
al. (1998). 



A /A* = 1 by this temperature. Although this could account 
for the coronal radiation losses, it is significantly smaller 
than the estimated energy loss by thermal conduction at 
1.4 x 10 6 K (~ 3 x 10 6 erg cm~ 2 s -1 , according to Jordan et 
al. 2001a). Also, there is insufficient flux below ~ 2 x 10 4 K. 
Thus the simple acoustic model appears to be unable to 
provide enough energy to heat the corona of t Eri. 

For an Alfven wave flux the argument concerning the 
area factor fails because the factor b in equation (12) can 
only be constant or decrease with T e . There is ample energy 
flux to account for the energy losses from the corona if b — 1 
(i.e. all the surface magnetic flux extends to the corona). 
Indeed, only around 20 per cent of the surface magnetic 
flux need extend to the corona to match the estimated en- 
ergy flux lost by thermal conduction at, and radiation losses 
above, 1.4 x 10 6 K (~ 4 x 10 6 erg cm~ 2 s -1 ). With the 
fluxes shown in Fig. 5 difficulties occur with the decrease in 
FntA/A* above 10 5 K. This decrease could be interpreted 
as being due to energy dissipation in the upper transition 
region, which could in principle be balanced by radiation 
losses and thermal conduction. However, including a heat- 
ing process in addition to radiation losses will act to increase 
the gradient of the theoretical emission measure distribu- 
tion (EMD) (see Jordan 2000). To match the decrease in 
F^tA/A*, apparent above 10 5 K in Fig. 5, would require 
the gradient of the observed EMD to be at least an order of 
magnitude greater than that derived from line fluxes. Thus 
an interpretation of the fluxes in Fig. 5 in terms of heating 
in the upper transition region is not viable. If b is less than 
1, then this constraint is less demanding. 

The radiation losses shown in Figs. 5 and 6 ignore the 
role played by thermal conduction above T e ~ 2 x 10 5 K. 
If the thermal conduction leaving the corona is balanced by 
the local radiation losses, then no additional heating is re- 
quired between 2 x 10° K and the base of the corona. Over 
this temperature range, when comparing with the heating 
flux, one needs to consider only the conductive flux and ra- 
diation lost from the corona, with further radiation losses 
being added below 2 x 10 5 K. 

The heating requirements of the chromosphere can also 
be considered. First we consider an acoustic flux. The form 
of Fig. 6 shows that if the turbulent velocities are associated 



with propagating waves a substantial fraction of the total 
non-thermal energy would be deposited in the low chromo- 
sphere. The non-thermal energy flux would have to decrease 
by ~ 10 9 ergs cm~ 2 s _1 across the low chromosphere. This 
far exceeds any reasonable estimate of the radiative losses 
from this region. Linsky & Ayres (1978) estimated that in 
cool main-sequence stars, the radiative fluxes in the Mg II h 
& k lines are comparable to 30 per cent of the total radiative 
flux from the chromosphere. The total observed flux at the 
Earth in the Mg II h and k lines in e Eri (in our STIS spec- 
tra) is around 5.4 x 10~ n erg cm~ 2 s _1 , corresponding to a 
stellar surface flux of 1.7 x 10 6 erg cm~ 2 s _1 . This gives an 
estimated total surface radiative flux from the chromosphere 
of 8.8 x 10 ergs cm - s _1 . Although only a crude estimate, 
this shows that far too much energy would be dissipated in 
the low chromosphere if the observed non-thermal motions 
are interpreted in terms of a propagating acoustic wave flux. 

If we invoke an area filling factor, then in the chromo- 
sphere, Aj A t . would have to be much lower than at 10 4 K, 
and would need to increase rapidly by around 10 4 K. This 
seems unlikely, as in the Sun, acoustic waves are thought to 
be important only in the interiors of supergranulation cells, 
whose area does not change significantly in the overlying re- 
gions. By analogy with the solar chromosphere, the spatially 
integrated emission and the observed line widths are likely to 
be dominated by the contribution from the supergranulation 
boundaries. In this case, relating the line widths to acous- 
tic waves will not be correct, and the apparent discrepancy 
between the observed and deduced fluxes is understandable. 

From the limits to the Alfven wave flux shown in Fig. 5 
it is clear that there is no difficulty in principle in obtaining 
sufficient flux to pass through to the lower transition region 
and the flux around 2 x 10 4 K represents an upper limit avail- 
able to the overlying regions. We do not address the question 
of how the supergranulation boundary regions in the chro- 
mosphere are heated. (Osterbrock 1961 gives a discussion of 
the possible processes involved, but with numerical values 
of parameters that require revision in the light of modern 
measurements of magnetic fields and current chromospheric 
models) . 

In Fig. 7 we show the non-thermal velocities calculated 
using the solar line widths given by Chae et al. (1998). Fig. 
8 shows the non-thermal wave fluxes calculated from equa- 
tions (12) and (14) using these non-thermal velocities and 
densities taken from model C (average quiet Sun) of Ver- 
nazza, Avrett, & Loeser (VAL) (1981). Again, the upper and 
lower limits on the Alfven wave fluxes are shown. As in Fig. 
5, b has been set to 1.0. The adopted temperatures of line 
formation are similar to those used by Chae et al. (1998), ex- 
cept for the neutral chromospheric lines. For these lines sig- 
nificantly lower values of Tform have been used, based on ra- 
diative transfer calculations with the VAL-C (1981) model. 
These lower temperatures lead to a difference in the depen- 
dence of the fluxes with T e , when comparisons are made 
with the calculations by Chae et al. (1998) (their fig. 13). 
The mean magnetic flux density was taken as 20 Gauss (see 
discussion in Montesinos & Jordan 1993). Apart from the ab- 
solute scales, the variations of the spatially averaged fluxes 
are quite similar to those found for e Eri: for Alfven waves 
the lower limit to the chromospheric flux is similar to the 
flux in the lower transition region, while for acoustic waves 
and the upper limit to the Alfven wave flux, there is a large 
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drop in Fnt between the chromosphere and transition re- 
gion. In both cases the flux is almost constant through the 
low transition region with a decrease in the upper transition 
region. In the Sun, the factor by which the flux decreases in 
the upper transition region is larger than in e Eri, but be- 
gins at a slightly higher temperature. The presence of a dip 
at around logT e — 4.5 relies heavily on the accuracy of the 
non-thermal contribution to the line-width of the optically 
thick C ii lines deduced by Chae et al. (1998). 

It is surprising that in the Sun, and in e Eri, in the 
low transition region £ increases according to T}J A , consis- 
tent with a constant Alfven wave energy flux in the WKB 
approximation. This approximation is valid only if the wave- 
length is less than the scale-height of the Alfven wave speed 
(see e.g. Meyer 1976; Hollweg 1991 and references therein); 
i.e. if the Alfven wave period t w is such that 

t w « — 2H (15) 

V a 

where H is the scale-height of the Alfven speed. If t w is much 
larger than this value, then the waves become evanescent 
and £ would be constant with height (which is observed to 
occur around 10 5 K). Defining H — v a (dh/dv a ), the critical 
period can be expressed as 



2.0 x l(T 10 £;m(0.3)iT e 3/2 



where the temperature gradient has been written in terms 
of the true (intrinsic) emission measure and the electron 
pressure (P e = N e T e ), i.e. 

PeP H dh 



dlogB 

5 T e ■ 



(16) 



Em(0.3) t 



(17) 



V2T e dT e ' 

For a plane parallel layer Em(0.3) t can be replaced by 
2Em(0.3)obsA*/A. As before, BA is replaced by B A*b, and 
the unknown area A, but not its variation with T e , cancels. 
Working above ~ 2 x 10 4 K (since we take P e and \i to be 
constant) the critical period can be found as a function of 
T e , provided assumptions are made about b and the change 
in the area factor A. I.e. 

2.0 x 10" 10 Sm(0.3) oi , s r ( ? /2 



B o 6P e 3/2 [0.5 + ™^] 



dlog(6/A) i 



(18) 



Since the area occupied by the solar supergranulation net- 
work varies little through the low to mid transition region, 
we assume that A is constant and we have to assume that 
6=1. For the Sun, using B = 20 G, P e = 2.1 xlO 14 cm -3 K, 
and the network values of Em(0.3) o bs from Macpherson & 
Jordan (1999), the values of r c are ~ 20 s at 2 x 10 4 K, 
~ 4 s at 8 x 10 4 K, and increase again to ~ 30 s at 3 x 10 5 K. 
For the assumed parameters these are minimum values since 
only a small increase in A with T e will increase r c , and if 
6 is less than 1, and decreases with T e , this will also in- 
crease t c . Also, the value of Bo used could be a factor of 2 
smaller. Thus although the critical wave periods are short, 
it is not certain that WKB waves can be excluded. The most 
stringent condition occurs around 8 x 10 4 K, where (maybe 
coincidentally) £ begins to increase by a power of T e which 
is less than 1/4. In e Eri, the values of Em(0.3) o bs, Pe and 
Bo are all larger, but result in values of r c that are smaller 
by a factor of about 130, with the same trends with T e . The 
values of b and dlog(6/ J 4)/dlogT e , which act to increase 



the critical period, therefore become the crucial factors in 
deciding whether or not the waves can satisfy the WKB ap- 
proximation. 

As Hollweg (1991, and references therein) has stressed, 
if the WKB approximation does not hold, then the energy 
flux carried can be larger or smaller than that deduced by us- 
ing the approximation. There is also the possibility of wave 
reflection by the steep Alfven wave gradient, and wave prop- 
agation at particular frequencies (resonances). Together, in 
principle, these effects might account for the apparent de- 
crease in the Alfven wave flux that appears in Figs. 5 and 
8, without invoking wave dissipation in the upper transition 
region. If wave reflection occurs, there is also the possibil- 
ity that eventual wave dissipation could contribute to the 
heating of the lower transition region. Although we do not 
attempt an analysis of possible wave processes, the new ob- 
servations of e Eri, and the new models being produced, will 
place strong constraints on future detailed calculations of 
wave propagation through the transition region. 



4 LINE SHIFTS 

Unless otherwise stated, the rest wavelengths given in Table 
3 are taken from Kurucz & Bell (1995). In general, these 
agree with those given by Kelly (1987) to within 0.001 A. 
In Jordan et al. (2001a) we noted that the wavelength in- 
terval between the two lines of Nv is slightly larger than 
that given by the values from Kelly (1987) (which are the 
same as in Kurucz & Bell 1995, and originate from labora- 
tory measurements by Hallin 1966). There has been previ- 
ous debate concerning the Nv wavelengths in the Sun (see 
Achour et al. 1995; Brekke, Hassler & Wilhelm 1997, and 
references therein). Brekke et al. (1997) find that for the 
1242. 8-A line the best agreement is with the measurement 
by Hallin (1966) (1242.804 A). For the 1238.8-A line, they 
find best agreement with the earlier measurement by Edlen 
(1934) (1238.800 A). Combining these wavelengths gives a 
wavelength interval of 4.004 A. Of the various combinations 
of the laboratory wavelengths this one does give the best 
agreement with our interval of 3.994 A, which should be 
accurate to within 0.007 A, and appears to be the best mea- 
surement of the interval to date. Since we do not have ab- 
solute wavelengths, we use the values by Hallin (1966) in 
Table 3, but point out that the 1242.804-A line is likely to 
give more accurate wavelength shifts. Our measured inter- 
val suggests that the wavelength of the other line should 
be 1238.810 A. It would be useful to remeasure these wave- 
lengths in the laboratory with modern techniques. 

New interferometric measurements of lines from the 
alkali-like ions C IV and Si IV are available from the work of 
Griesmann and Kling (2000) , and these are adopted in Table 
3. For C IV the new values agree more closely with the values 
of Kelly (1987) than with those in Kurucz & Bell (1995) and 
lead to smaller observed red-shifts. The wavelength inter- 
val between the two C IV lines that we measure is definitely 
smaller than that given in Kurucz & Bell (1995), and agrees 
most closely with that from the wavelengths given by Kelly 
(1987). For the Siiv lines our measured interval is consistent 
with those listed in all the above three papers. The new mea- 
surement by Griesmann & Kling (2000) for the line of Alii 
is also adopted in Table 3, although it is not significantly 
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Figure 8. logFNT versus logTf orm from the solar data of Chae et. al (1998). The left-hand panel is for Alfven waves and the right-hand 
panel is for acoustic waves. In the left-hand panel the open symbols indicate the calculations using pi on and the filled symbols show 
those with the total density p. 



different from the earlier value, given the uncertainty in the 
observed value. 

The wavelength calibration of FUSE is not accurate, so 
no absolute wavelengths can be measured for the O VI and 
C n 1037-A mult, lines. Since these lines are all close together 
in wavelength, their relative wavelengths can be measured. 
The wavelength interval between the O VI lines agrees with 
the values in both Kurucz & Bell (1995) and Kelly (1987) 
to within our measurement uncertainty of 0.025A. That be- 
tween the C II lines agrees with the laboratory measurements 
for this spectroscopic standard, to within our measurement 
uncertainty. There is then no relative shift between the O VI 
lines and the C II lines (which are formed where the mean 
red-shift is 2-7 km s _1 ), to within the measurement uncer- 
tainty of ~ 4 km s~\ it appears that the red-shifts of the 
O VI lines are no greater that that of O v. 

The line shifts given in Table 3 are plotted as a func- 
tion of wavelength in Fig. 9 and of temperature of formation 
in Fig. 10. The diamonds indicate lines observed with the 
E140M grating and the asterisks those which are observed 
with the lower resolution E230M grating. For multiplets in 
which the individual line shifts agree to within the measure- 
ment errors (see Table 3) the individual shifts have been 
combined to give a single shift for the multiplet. 

Fig. 9 shows that most of the lines display a small 
(few km s _1 ) red-shift, although the scatter is quite large. 
There appears to be an offset in the wavelength calibration 
of the E230M grating (which covers the wavelength range 
A > 1750 A) compared with that of the E140M grating, since 
at a given T e , the lines observed around 1800 A to 1900 A ob- 
served with the E230M grating show larger redshifts than 
those observed with the E140M grating. The offset is small 
and corresponds to only about 2 pixels in wavelength, but in 
the discussion below we restrict our comments to the trends 
in the E140 data. 

Fig. 10 does not show a strong trend in the line shifts 
with Tform- Using only the lines observed with the E140M 
grating, the mean red-shift is 3.7±0.6 km s _1 . The mean red- 
shift of the mid transition region lines (4.4 < log Tf orm < 5.5) 
is 4.9±0.8 km s _1 while that for the lower temperature lines 
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Figure 9. Line shift versus wavelength. The lines observed with 
the E140M grating are shown by o, and those observed with 
the E230M grating by *. 
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Figure 10. Line shift versus logTf orm . Symbols as in Fig. 9. 
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Table 4. Fit parameters from two-Gaussian fits to selected strong lines. F denotes the integrated line flux, given in 10 -13 ergs cm~ 2 s _1 
at the Earth. The subscripts N and B denote "Narrow" and "Broad" components respectively. The errors are formal la errors deduced 
from the correlation matrix of each fit. 



Ion 


Ao(A) 


F-n 


AA N (A) 


^obs,N(A) 


Fb 


AAb(A) 


A bs,B(A) 


C iv 


1548.204 


3.293 ±0.152 


0.165 ±0.004 


1548.221 ± 0.001 


2.177 ±0.137 


0.420 ±0.019 


1548.247 ± 0.004 


C iv 


1550.781 


1.396 ±0.108 


0.147 ±0.007 


1550.784 ± 0.002 


1.453 ±0.103 


0.372 ± 0.016 


1550.805 ± 0.004 


N v 


1238.821 


0.676 ± 0.062 


0.151 ± 0.008 


1238.828 ± 0.002 


0.304 ± 0.058 


0.377 ± 0.047 


1238.845 ± 0.009 


N v 


1242.804 


0.334 ±0.051 


0.147 ±0.011 


1242.817 ±0.003 


0.145 ±0.047 


0.332 ± 0.062 


1242.865 ± 0.028 


O vi 


1031.912 


3.540 ±0.136 


0.173 ± 0.003 


1031.940 ±0.001 


1.065 ±0.127 


0.376 ± 0.022 


1031.959 ±0.005 


O vi 


1037.613 


1.858 ±0.074 


0.170 ± 0.004 


1037.615 ±0.001 


0.417 ±0.067 


0.392 ± 0.035 


1037.630 ± 0.006 


Si in 


1892.030 


0.794 ± 0.089 


0.158 ± 0.012 


1892.105 ±0.004 


0.523 ±0.093 


0.427 ± 0.042 


1892.112 ±0.022 


Si iv 


1393.760 


1.285 ±0.050 


0.136 ± 0.004 


1393.780 ± 0.001 


0.960 ± 0.047 


0.379 ± 0.012 


1393.793 ± 0.003 


Si iv 


1402.773 


0.739 ± 0.038 


0.141 ± 0.005 


1402.794 ±0.002 


0.450 ± 0.034 


0.427 ± 0.026 


1402.803 ± 0.006 



is 2.1± 1.1 km s _1 . Thus, on average, the lower temperature 
lines are less red-shifted by about 3 km s _1 . The Fe XII 
line is apparently blueshifted, but this relies on the solar 
wavelength at the limb. 

In the solar transition region, using spatially integrated 
observations at Sun-centre, Peter & Judge (1999) found red- 
shifts of up to 9.9 km s~ x in the Nv lines, and smaller 
red-shifts in lines formed at lower and higher temperatures. 
These red-shifts decrease towards the limb. For the Siiv, 
C iv and N v lines in common with our observations of e Eri, 
we find similar red-shifts (if we adopt the wavelengths from 
Kurucz & Bell 1995), except for the Nv lines, for which 
our values are significantly smaller. The new wavelengths 
we adopt for Siiv and C iv do, however, lead to smaller red- 
shifts (see Table 3) . If the Fe XII wavelength is reliable, then 
like Peter & Judge (1999) we observe a blue-shift in the 
Fexn line, although this is smaller in e Eri. 

The physical origin of the shifts is not yet understood, 
although it has been suggested by Hansteen (1993) that 
these could be due to waves propagating downwards from 
reconnection events (nanofiares) in the corona, or in the mid 
transition region (Peter & Judge 1999). However, our line- 
width measurements, using one-Gaussian fits, show no ev- 
idence that the energy flux increases with temperature, as 
one would expect if the origin of the waves was within the 
corona. In the solar photosphere the velocity fields associ- 
ated with the supergranulation flow are very small, but the 
possibility of convective overshoot into the transition region, 
and resulting downflows in the cell boundaries, does not ap- 
pear to have been investigated. With lower densities and 
mass conservation, larger velocities are in principle possible 
at the level of the transition region. It is also interesting 
that, in at least one solar boundary region, Macpherson, 
Jordan & Smith (1999) found that the downflows appear to 
be adjacent to, rather than in, the brightest regions of the 
boundaries. Overall, the origin of the red-shifts is still not 
understood. 

As discussed above, in stellar observations only the com- 
bined effects of red- and blue-shifts can be observed, which 
could contribute to the line widths. The shifts observed are 
significant compared with the uncertainty of 0.9 km s _1 
in the stellar radial velocity. The shifts derived from two- 
component fits to line profiles are discussed below in Section 
5. 



5 TWO-COMPONENT FITS 

It has previously been found (see Section 1) that the pro- 
files of some transition region lines are described better by 
a two-Gaussian fit than a one-Gaussian fit. This has been 
taken as an indication that there are (at least) two separate 
atmospheric components contributing to the line emission. 
This two-component nature of line widths has been stud- 
ied in the Sun (Peter 2000a,b; 2001) and other stars (Wood 
ct al. 1997). These authors have found that line profiles for 
strong lines are described well by a sum of two Gaussians, 
one narrow and one broad. For the Sun, Peter (2000a,b; 
2001) found that the presence of the second, broad com- 
ponent is restricted to the bright supergranulation network 
boundaries and that the cell interior regions show only the 
narrow component. 

The excellent resolution and high signal-to-noise of the 
e Eri STIS spectra allow reliable two-component fits to be 
made to several of the strong spectral lines. The lines for 
which two-component fits have been made (these were se- 
lected on the grounds that they are sufficiently strong to 
reliably identify a second component) are listed in Table 4, 
which also gives the flux (F), wavelength (A b s ) and width 
(AA) for both the narrow and broad components. The er- 
ror estimates quoted in Table 4 are formal la errors derived 
from the covariance matrix of the fit. The fluxes and widths 
are all strongly correlated while (with the exception of the 
N v 1242-A line) the wavelengths are only weakly correlated 
with the other parameters. Errors estimated this way do not 
explicitly include the limitation on measuring wavelengths 
due to the instrumental resolution, so the true uncertain- 
ties in A obs are larger than estimated here for the narrow 
components, and are similar to the errors in A bs as given 
in Table 1. Based on the line lists of Kurucz & Bell (1995) 
and Sandlin et al. (1986), no significant blends are expected 
in the wavelength ranges over which the fits to the profiles 
have been made. 

For the strong lines, a two-component fit gives a sig- 
nificantly better description of the data than does a one- 
component fit. To illustrate this, Fig. 11 shows the ob- 
served Civ, Siiv, Nv, Ovi and Sim (1892- A) lines with 
single-Gaussian fits and residuals (left-hand panels) and 
with double- Gaussian fits and residuals (right-hand pan- 
els). By examining the residuals it is clear that the two- 
component fit gives a substantial improvement for the C IV 
and Si IV lines and some improvement for the other lines. To 



Velocity fields in e Eri 13 





Figure 11. Comparison of one-Gaussian fit (left) and two-Gaussian fit (right) for (from top left) the Civ 1548-A and 1550- A lines, 
the Si IV 1393-A and 1402-A lines, the Nv 1238- A and 1242-A lines, the Ovi 1032- A and 1038-A lines and the Sim 1892-A line. The 
observed profiles are shown by the solid line, the fits by the dashed lines and the residuals by the dotted lines. 
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Table 5. Reduced chi-squared (x r ) f° r one " an< l two-Gaussian 
fits to selected strong lines. 



Ion 


Ao (A) 


Xr (one-Gaussian) 


Xr (two-Gaussian) 


C iv 


1548.204 


14.28 


0.87 


C iv 


1550.781 


8.24 


1.40 


N v 


1238.821 


2.18 


0.80 


N v 


1242.804 


2.20 


1.22 


O vi 


1031.912 


5.25 


1.34 


O vi 


1037.613 


3.11 


1.18 


Si in 


1892.030 


1.70 


0.89 


Si iv 


1393.760 


18.04 


2.18 


Si iv 


1402.773 


6.67 


1.06 



be more quantitative, the one- and two-Gaussian fits have 
been used to calculate the reduced chi-squared {xV) statistic 
for each of the lines considered; these are given in Table 5 
and show that the two-Gaussian fits give a significantly bet- 
ter description of the profiles than do the one-Gaussian fits. 
The reduction in \r is greatest for the C iv and Si IV lines 
but is still large enough to strongly suggest the presence of 
a second component in the Nv, O vi and Sim lines. 

As for the one-Gaussian fit parameters, the two- 
Gaussian fit parameters (and associated la errors) have 
been converted to velocity shifts (v) and most probable tur- 
bulent velocities (£) for both components of each line. The 
results are given in Table 6. The temperatures at which both 
components are formed have been assumed to be the same 
as that used for the one-Gaussian fits. The uncertainty in £ 
allows for a ±0.1 dex uncertainty in logTf orm in the same 
way as adopted for the one-Gaussian fits. For O vi no ve- 
locity shift is given for the narrow component owing to the 
uncertainty in the FUSE wavelength calibration. The veloc- 
ity shift given for the broad component is that relative to 
the narrow component. 

The red-shifts found for the narrow components are es- 
sentially the same as those found for the one-Gaussian fits. 
A significant difference between the present results and the 
solar results found by Peter (2000a,b; 2001) is that in e Eri 
the broad component is red-shifted relative to the narrow 
component, typically by about 4 km s _1 . This can be seen 
directly in the asymmetry of the residuals shown in Fig. 11. 
Peter (2000a) found that in observations of C iv near Sun- 
centre the narrow component in the supergranulation cell 
boundaries is usually red-shifted, but the broad component 
is red-shifted less. However, in Peter's (2000a) detailed plots, 
there are a significant number of points where the red-shift 
of the broad component exceeds that of the narrow com- 
ponent, although these tend to occur when the intensity of 
the broad component is small. Since the relative intensity of 
the broad component in e Eri is larger than in the Sun (see 
below), this trend acts in the wrong direction to account 
for our results. In the cell interiors, where the broad com- 
ponent is absent, the profiles show red-shifts comparable to 
those of the narrow components in the boundaries. In stellar 
observations of the Civ lines Wood et al. (1997) observe a 
relatively greater red-shift of the broad component in the 
M-type flare star AU Mic, but this line was not observed in 
the other less active main-sequence stars a Cen A and B. 
However, in all three stars the Siiv lines show relative shifts 
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Figure 12. log£ versus logTf orm for the 2-componcnt fits (broad 
components shown by triangles; narrow components shown by di- 
amonds). The single-component fits are shown by squares. The 
full line shows the sound speed and the dashed line shows a gra- 
dient of £ oc TfJ^ , with an amplitude normalized to the O VI 
narrow component. 

in the same sense as in the Sun, in contrast to our results 
for e Eri. 

Fig. 12 shows the non-thermal velocities of the two com- 
ponents as a function of the line formation temperature. For 
comparison, the single component points are shown again. 
The solid line shows the sound speed. The dashed line shows 
a gradient corresponding to £ oc T f *^, normalized to pass 
through the O vi narrow components. 

The width of the broad component gives non-thermal 
velocities which are greater than the sound speed for the 
Si in and Si iv lines, close to the sound speed for C iv and 
significantly below the sound speed for N v. This behaviour 
is similar to that found by Peter (2001) for the Sun and for 
the high gravity stars studied by Wood et al. (1997). The 
widths of the low temperature lines from the single-Gaussian 
fits are also close to the sound speed (see Fig. 1 for the error 
bars) . 

As in the Sun (Peter 2001), the ratio of the non-thermal 
widths of the two components is almost constant. The mean 
ratio in e Eri is about 3, larger than the solar value (for lines 
in common) of about 2. Although the non-thermal widths 
of the broad components in e Eri are on average similar to 
those in the Sun, the widths of the narrow component are 
smaller in e Eri, by about a factor of 1.3. Wood et al. (1997) 
found a systematic decrease in the width of the C iv and Si iv 
narrow components with increasing surface gravity. These 
widths in e Eri fit their proposed relation £ oc g~ 0S8±0 ' 0r 
remarkably well. A dependence on g t , or alternatively the 
gas pressure, could indicate that turbulent motions are ul- 
timately involved in any wave generation. Peter (2001) sug- 
gests that the broad component corresponds to a constant 
Alfven wave flux, with £ oc Tf/^. In e Eri, to within the 
uncertainties, the same form would fit both components up 
to T e ~ 3 x 10 5 K. As stressed above, the temperature de- 
pendence expected for £ depends on the details of the wave 
periods and any departures from the WKB approximation. 

Peter (2001) found that the relative fluxes in the nar- 
row and broad components vary with temperature, with 
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Table 6. Line shifts (v) and turbulent velocities (£) for two-component fits. The subscripts N and B denote "Narrow" and "Broad" 
components respectively. The errors are deduced from the Icr errors in Table 4 and are given in same format as those in Table 2. 



Ion Aq (A) i>n (km s - 1 ) £n (km s - 1 ) db (km s - 1 ) £b (km s 1 ) 



C iv 


1548 


204 


+3.3 ±0.2 


15 


ql +1.23(0 
-1.52(0 


69) 
92) 


+8.3 + 0.8 


47 


7 „+2.49(0.24) 
-2.58(0.30) 


C iv 


1550 


781 


+0.6 ±0.4 


13 


97 +1.80(0 
Z '-2.29(l 


82) 
11) 


+4.6 ±0.8 


41 


oo+2.17(0.27) 
-2.28(0.34) 


N v 


1238 


821 


+1.7 ±0.5 


16 


q „+2.50(l 
-3.19(1 


10) 
49) 


+5.8 ±2.2 


53 


nc . + 7.36(0.36) 
U -7.64(0.46) 


N v 


1242 


804 


+3.2 ±0.7 


16 


,0+3.09(1 
id -4.04(l 


15) 
58) 


+14.7 ±6.8 


46 


nq +9. 67(0.42) 
u -10.16(0.53) 


O vi 


1031 


912 


see text 


18 


„ 9 +2.30(l 
-3.09(2 


53) 
12) 


+5.5 ± 1.7 


61 


„„+4.54(0.49) 
-4.80(0.62) 


O vi 


1037 


613 


see text 


17 


7 „+2.66(l 
-3.65(2 


62) 
28) 


+4.3 + 2.0 


63 


„ fi +6.87(0.47) 
-7.19(0.59) 


Si in 


1892 


030 


+11.9 ± 0.6 


13 


,. + 1.40(0 
-1.50(0 


13) 
16) 


+13.0 + 3.5 


40 


n .+4.10(0.04) 
U -4.13(0.05) 


Si iv 


1393 


760 


+4.3 ±0.2 


16 


„ n +0.78(0 
-0.87(0 


23) 
30) 


+7.1 + 0.6 


48 


, 7 +1.64(0.08) 
-1.67(0.10) 


Si iv 


1402 


773 


+4.5 ±0.4 


16 


„ 7 +0.91(0 
-0.99(0 


23) 
29) 


+6.4 + 1.3 


54 


.-+3.43(0.07) 
-3.46(0.09) 



Table 7. The ratio of the fluxes in the two-components including 
the lcr error estimates. The subscripts N and B denote "Narrow" 
and "Broad" components respectively. 



Ion 


Ao (A) 




C iv 


1548.204 


0.67 + 0.07 


C iv 


1550.781 


1.04 + 0.15 


N v 


1238.821 


0.45 + 0.13 


N v 


1242.804 


0.43 + 0.21 


O vi 


1031.912 


0.30 + 0.04 


O vi 


1037.613 


0.22 + 0.05 


Si in 


1892.030 


0.66 + 0.13 


Si iv 


1393.760 


0.75 ± 0.06 


Si iv 


1402.773 


0.61 ± 0.07 
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the broad component having its largest contribution at 
logTform w 5.1 with smaller contributions at both lower and 
higher temperatures (~ 0.2 for Siiv and Ovi and ~ 0.4 for 
C iv). The ratio of the fluxes in the two components in e Eri 
are given in Table 7. The la errors on the flux ratio are 
somewhat larger than might be expected given the errors in 
Table 4. This is because of the strong anti-correlation be- 
tween Fb and Fs which has been accounted for in the la- 
error estimates given in Table 7. In e Eri the flux in broad 
component is significantly larger, relative to the narrow one, 
than in the Sun. If the broad and narrow components gen- 
uinely represent different physical regions, this suggests that 
the broad component may increase with stellar activity. This 
was proposed by Wood et al. (1997) on the basis of a larger 
sample of stars for which a correlation with the stellar X-ray 
flux was apparent. 

Fig. 13 shows a plot of the ratio of the line fluxes in 
the two components versus logTform- Although there is a 
significant difference between the flux ratio deduced from 
the two Civ lines, the Nv and Ovi lines do show that the 
broad component is less important at higher temperatures. 
In contrast to the results of Peter (2001, his fig. 6), in e Eri 
the Si iv and Sim lines do not show a rapid decrease in 
the contribution of the broad component at lower temper- 



Figure 13. Ratio of the flux emitted in the broad component to 
the that in the narrow component versus log Tf orm . 



atures. We point out that the sharp drop in the relative 
strength of the broad component at low temperatures de- 
duced by Peter (2001) depends on his analysis of the Sill 
1533-A and C n 1335-A lines (the latter is actually a close 
blend of two components of the same multiplet). Both of 
these lines have been excluded from the current study on 
the grounds that they have significant opacity and hence 
their profiles can be affected by radiative transfer effects. 
Also, models of the atmosphere (Sim 2002) suggest that the 
temperatures of formation adopted by Peter (2001) for these 
lines are too high, in which case the thermal broadening has 
been overestimated. The C n 1335-A line in e Eri does show 
broad wings and the analysis of this line profile will be dis- 
cussed in later work on chromospheric modelling, including 
the radiative transfer. 

Comparing Fig. 12 and 13 suggests a possible reason 
for the flattening off of the non-thermal velocities with T e 
above 10 J K when single Gaussian fits are used (see Fig. 1). 
Above this temperature the contribution of the broad com- 
ponent to the total line flux decreases so that the combined 
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Figure 14. log Fnt ^4/^4* versus logTf orm for 2-component fits 
(broad components marked with a triangle and narrow compo- 
nents with a diamond). The fluxes are more constant than when 
a single-Gaussian fit to the line profiles is used. The boxes indi- 
cate results using single-Gaussian fits, computed using pi on . (See 
Fig. 5 for results using p.) The full line shows the flux required 
to account for the radiation losses above a given T e . 



profile becomes dominated by the narrow component. Thus 
the single-Gaussian fit produces a smaller line width, even 
though the widths of both the narrow and broad component 
are not decreasing with logTf orm . This explanation depends 
on the Fe XII line flux being dominated by the narrow com- 
ponent, while the mid transition region lines Si III, Si IV, 
C IV and N v contain significant contributions from both 
components. 

The non-thermal flux in each component for waves trav- 
elling at the Alfven speed, and assuming the same P e in both 
components, is shown in Fig. 14. Only calculations using pi on 
are shown since the lines for which two- Gaussian fits have 
been carried out all form at sufficiently high temperatures 
that pion ~ p. The drop in the flux derived for the narrow 
component (spatially averaged as before) between 10° K and 
1.4 x 10 6 K is now only a factor of 1.8, compared with a fac- 
tor of 4.5 when single-Gaussian fits are made. The fluxes are 
still large compared with those required and the wave prop- 
agation needs to be investigated. We have no way of telling 
whether or not the flux in the broad component continues 
to higher temperatures. If it decreases to the value derived 
from the Fexn line, then there are problems in disposing of 
the energy, unless b and P e are small throughout the broad 
component region. There would be similar problems regard- 
ing deposition of the energy if the broad components were 
restricted to regions where the magnetic flux does not con- 
tinue to the corona. As stressed above the broad component 
flux could continue without being observable in the Fexn 
line profile. In this case, the correlation between the degree 
of the contribution of the broad component to the line pro- 
files and the X-ray flux, found by Wood et al. (1997) is 
surprising, since the radiative losses should be dominated 
by the narrow component. It is possible that the X-ray flux 
and the broad components are related to different aspects 
of the magnetic field configuration, which themselves have 
systematic trends. Thus the origin of the broad component 



and its relation to coronal heating is still unclear. Energy 
balance models are required to explore this further. 

At lower temperatures, extrapolating the results of Pe- 
ter (2001) would suggest that the lines are dominated by 
the narrow component. Certainly, there is no strong evi- 
dence for the presence of a second component in the C I 
or O I intersystem lines. If these lines were dominated by 
the broad component, the large rise in £ (see Fig. 12) when 
combined with pi on , would at face value imply an increase 
in the wave flux and an energy source in the low transi- 
tion region/upper chromosphere. For example, Axford & 
McKenzie (1992) have suggested that reconnection events 
involving small flux tubes at the level of the chromosphere 
or low transition region could generate Alfven waves. These 
could be high frequency waves which could propagate to 
the corona before being damped (McKenzie, Axford & Ba- 
naszkicwicz 1995, McKenzie, Banaszkiewicz & Axford 1997). 
If p is adopted, the wave flux in the chromosphere exceeds 
that required at higher temperatures (see Section 3) and it 
is possible that the broadening of the chromospheric lines is 
the result of a different heating process and not a continua- 
tion of either transition region component. 

An alternative explanation of the line-widths (of one or 
both components) is that the ion and electron temperature 
are not equal and the widths reflect high ion temperatures 
(see Peter 2001). The idea is attractive given observations 
of high O VI ion temperatures in polar coronal holes (see 
Cranmer et al. 1999 and references therein), and proposals 
by Marsch, Goertz & Richter (1982) and McKenzie et al. 
(1995, 1997) that damping of Alfven waves by ion-cyclotron 
resonance could heat and drive the solar wind. However, to 
predict the resulting effects on ions in the transition region 
requires detailed calculations. 



6 CONCLUSIONS 

The excellent properties of the STIS instrument have al- 
lowed us to determine the non-thermal motions and line 
shifts in the outer atmosphere of e Eri through analyses of 
the profiles of optically thin lines. The non-thermal velocity 
fields are important in semi-empirical modelling of the at- 
mosphere and the results presented here have been used in 
such models (Sim 2002). 

The measured line widths have been used to derive the 
non-thermal energy content of the atmosphere. The energy 
which could be carried in acoustic or Alfven waves has been 
investigated using the simple WKB approximation. The crit- 
ical wave periods for this approximation to hold have been 
expressed in terms of observable parameters, as far as is pos- 
sible. The unknown parameters are related to the magnetic 
field. With new models of the outer atmosphere more rigor- 
ous calculations could be made of how such waves may or 
may not propagate to the corona. As concluded from pre- 
vious studies, the observations appear to rule out acoustic 
waves as the source of coronal heating. 

A detailed analysis of the two-component nature of the 
mid-transition region lines of Sim, Siiv, Civ, Nv and Ovi 
has been carried out, and the results compared with those 
found for the Sun and other main-sequence stars. It is found 
that the broad component contributes a larger fraction of 
the transition region line flux than in the Sun, supporting 
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the suggestion by Wood et al. (1997) that this component 
to the flux increases with stellar activity. In agreement with 
analyses of solar data the broad component is found to be 
less important in the higher temperature lines of Nv and 
Ovi; it might not be significant in the Fexn line. If so, 
the correlation between broad transition region lines and 
stellar X-ray emission found by Wood et al. (1997) is hard 
to understand, unless both independently reflect systematic 
changes in the structure and strength of the magnetic fields. 
In e Eri, the broad component is red-shifted relative to the 
narrow component, whereas in the Sun, relative blue shifts 
are found. 

The narrow components are less broad than in the Sun 
and follow the correlation between £ and g* found by Wood 
et al. (1997). We consider that there are fewer difficulties 
in relating the narrow component to the passage of Alfven 
waves which could heat the corona, in contrast to Peter's 
(2001) proposal that this role is played by the broad com- 
ponent. 

A large amount of work remains to be done on the theo- 
retical interpretation of the observations. However, the qual- 
ity of the results from STIS places new constraints on pro- 
posed heating mechanisms. 
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APPENDIX: BROAD WINGS DUE TO 
NON-ISOTROPIC VELOCITY FIELDS 

The purpose of this appendix is to draw attention to the 
possibility that the broad wings in observed transition re- 
gion line profiles (Section 5) might be due to a non-isotropic 
non-thermal velocity field, rather than being the result of 
emission in two physically distinct atmospheric components. 

In deriving equation (1), it is assumed that the non- 
thermal velocities are described by a Maxwellian distribu- 
tion in three-dimensions. There is no direct observational 
evidence to justify this assumption and it is interesting to 
examine the consequences of it being incorrect. 

Suppose that the non-thermal motions are the result of 
oscillations caused by transverse waves propagating along 
magnetic field lines. Consider a small volume element of 
material which is threaded by magnetic field lines which 
are at an angle 9 to the observers line-of-sight. The non- 
thermal motions will be restricted to the plane perpendicu- 
lar to the magnetic field. Assuming, for simplicity, that the 
non-thermal motions within this plane are described by a 2- 
D Maxwell distribution, the observed profile of an emission 
line emitted by the material in the volume element will be 
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where A is the wavelength measured relative to line centre, 
7(A) is the intensity at wavelength A and C is a constant, a 
and b specify the widths of the 2-D Maxwell distribution of 
non-thermal velocities and the 3-D Maxwell distribution of 
thermal velocities respectively. Line broadening from sources 
other than thermal and non-thermal motions have been ne- 
glected. 

In order to obtain the profile that would be observed 
from the entire star it is necessary to assume a distribution 
for the angle 9. Here we adopt the simplest assumption, 
that the distribution is spherically symmetric. This would 
be a reasonable assumption for a star if the magnetic field 
were completely disordered (randomly oriented) on macro- 
scopic scales. It would also be approximately valid if the field 
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Figure Al. Profile computed from Equation (A2) (solid line) 
with Gaussian over-plotted (broken line) and residual (dotted 
line). The values a = 0.15 A and b = 0.05 A have been adopted. 



was purely radial over an (unresolved) hemispherical stellar 
surface: the calculations presented here are appropriate for 
emission from a hemisphere but they do not account for the 
excess emission off the limb from material in the back hemi- 
sphere of the star. If emission from the back hemisphere were 
included in the calculation it would enhance the strength of 
the broad wings predicted in the 2-D case. The spherically 
symmetric distribution will also be appropriate for spatially 
resolved observations of the Sun, provided that the field is 
tangled on scales less than the spatial resolution or contains 
both radial and closed loop field configurations. If the field 
were locally ordered, then it would be necessary to adopt 
a more complex distribution. Adding up different 9 angle 
strips from centre to limb, the total stellar profile in the 
presence of 2-D non-thermal motions is 
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Alternatively, consider the case of non-thermal motions 
generated by longitudinal waves propagating along the mag- 
netic field. The non-thermal velocities are then confined to 
only one dimension, along the field. By analogy with the 
case given above the total stellar profile is 

7(A)i_ D oc f Sm 6 de = exp(-A 2 /(b 2 + a 2 cos 2 9)) . 

Jo v o + a cos 2 9 

(A3) 



Figs. Al and A2 shows profiles calculated from equa- 
tions (A2) and (A3), respectively. In both cases the ther- 
mal width has been chosen to be appropriate to the C IV 
lines, b = 0.05 A. The non-thermal width parameter (a) 
was chosen arbitrarily (see below). In each figure a single- 
Gaussian fit to the line and the residual for that fit are plot- 
ted. These should be compared with the fits and residuals 
for the observed profiles shown in Fig. 11. Two-Gaussian fits 
have been performed on the profiles calculated from equa- 
tions (A2) and (A3). The results of these fits (fluxes and 
widths for the two components) are given in Table Al. 

The values of a used in the calculations were chosen 
so that the width of the narrow Gaussian fit component 
is comparable to that observed for C IV. For the 2-D case 
a = 0.15 A and for the 1-D case a — 0.32 A. 

Figs. Al and A2 show that these non-isotropic velocity 
distributions produce line profiles that have broader wings 
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Figure A2. Profile computed from Equation (A3) (solid line) 
with Gaussian over-plotted (broken line) and residual (dotted 
line). The values a = 0.32 A and b = 0.05 A have been adopted. 



Gaussian profiles whose properties are similar to those of 
observed transition region lines. 



Table Al. Two-Gaussian fit to the computed profiles for a 2-D 
non-thermal velocity distribution (equation A2) and a 1-D non- 
thermal velocity distribution (equation A3). The fluxes (F) are 
given in arbitrary units and line widths (AA, FWHM) are given 
in A. The subscripts N and B refer to the narrow and broad 
Gaussian components. 



F N AAjv F B A\ B 



2-D case 0.0571 0.137 0.1560 0.244 
1-D case 0.0933 0.132 0.1253 0.394 



than a Gaussian, in qualitative agreement with the observed 
profiles. The broad wings are relatively stronger in the 1-D 
than 2-D case. The width of the broad Gaussian component 
fitted to the calculated profile is in very good agreement with 
that fitted to the observed C iv profiles for the 1-D case, but 
is too narrow in the 2-D case. There is some quantitative dis- 
agreement between these calculations and the observations. 
In both cases, the flux ratio of the broad to narrow compo- 
nent is too large, compared with the observations. The 1-D 
flux ratio is consistent with the observed 1550-A flux ratio 
at the 2-<t level, but is very different from the 1548-A ratio. 
Also, the computed profiles are symmetric about line centre 
while the observed broad components are all redshifted rel- 
ative to the narrow components by several km s _1 . In the 
model used here (which neglects the line formation process) , 
to produce a redshifted broad component would require a 
mean net flow away from the observer. To shift the broad 
component relative to the narrow for the 1-D case would 
require a redshift that is largest in the regions with small 9, 
while for the 2-D case the redshift would need to be largest 
in regions with large 9. 

In conclusion, it is apparent from Figs. Al and A2 that 
a non-isotropic velocity distribution could predict line pro- 
files with broad wings without the need to invoke a physi- 
cally distinct second component to the emission. Although 
the discussion here has used the passage of transverse and 
longitudinal waves as examples, it is noted that the results 
are general to any situation where the non-thermal veloc- 
ities are confined to only one or two-dimensions, whatever 
their origin. It is not proposed that the extremely simple 
one and two-dimensional non-thermal velocity distributions 
used here are realistic; we wish only to draw attention to 
the fact that such distributions can naturally lead to non- 



